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We reconsider the strong-coupling expansion for the Hubbard model recently introduced by 
Sarker and Pairault et al. By introducing slave particles that act as projection operators onto 
the empty, singly occupied and doubly occupied atomic states, the perturbation theory around the 
atomic limit distinguishes between processes that do conserve or do not conserve the total number 
of doubly occupied sites. This allows for a systematic t/U expansion that does not break down at 
low temperature {t being the intersite hopping amplitude and U the local Coulomb repulsion) . The 
fermionic field becomes a two-component field, which reflects the presence of the two Hubbard bands. 
The single-particle propagator is naturally expressed as a function of a 2 x 2 matrix self-energy. 
Furthermore, by introducing a time- and space-fluctuating spin-quantization axis in the functional 
integral, we can expand around a "non-degenerate" ground-state where each singly occupied site 
has a well defined spin direction (which may fluctuate in time) . This formalism is used to derive the 
effective action of charge carriers in the lower Hubbard band to first order in t/U . We recover the 
action of the t-J model in the spin-hole coherent-state path integral. We also compare our results 
with those previously obtained by studying fiuctuations around the large-?7 Hartree-Fock saddle 
point. 

PACS Numbers: 71.10.Fd, 71.27.Fd, 71.30-)-h 
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I. INTRODUCTION 

In the last two decades, the discovery of heavy-fermion 
compounds, high- Tc superconductors and organic con- 
ductors has revived interest in strongly correlated elec- 
tron systems. Although the first investigations go back to 
the 1960s, a proper understanding of correlation effects in 
Fermi systems remains a fundamental issue in cpjidensed- 
matter physics. Even for the Hubbard model,ucl which 
is supposed to be one of the simplest (realistic) mod- 
els of strongly correlated fermions, exact solutions or 
well-controlled approximationSf-jexist only in a few spe- 
cial cases, like in one-dimensiorJj or in the limit of infinite 
dimension. tl 

In a system where the Coulomb interaction dominates, 
it is natural to treat the latter exactly and then consider 
the kinetic energy within perturbation theory. In the 
Hubbard model, this amounts to expanding around the 
atomic limit (no intersite hopping: i = 0). [This kind 
of expansion will be referred to as strong- coupling ex- 
pansion.] One of the first examples of a strong-coupling 
expansion is due to Anderson. Before the introduction of 
the Hubbard model, he showed that Fermi systems with 
a strong local Coulomb repulsion are described by an 
effective Heisenberg Hamiltonian at half-filling (one par- 
ticle per site on average). The strong on-site Coulomb 
repulsion U completely freezes out the kinetic energy, 
and virtual intersite hopping results in an effective an- 
tiferromagnetic (AF)-|ejx;change interaction with coupling 
constant J = 4t2/[/.tiy 

The strong-coupling expansion for the Hubbard model 



presents technical difficulties that do not show up in 
standard weak coupling perturbation theories. Since the 
atomic Hamiltonian is not quadratic (it contains the on- 
site Coulomb repulsion), there is no Wick's theorem so 
that the standard many-body techniques cannot be ap- 
plied. The perturbation theory in the hopping amplitude 
is based on a cumulant expansion with no linked cluster 
theorem. Given these technical difficulties, it is not sur- 
prising that contradicting results for the expansion of the 
free energy and other th^rpapdynamical quantities can be 
found in the liter ature.EThEl The calculation of dynami- 
cal quantities turns out to be even more involved. The 
first attempt to calculat£J:he single-particle Green's func- 
tion is due to HubbardED His approach became effective 
only with the latter dpjelopment of a Wick's theorem 
for Hubbard operators ila The strong-coupling expansion 
based on Hubbard operators remains however quite cum- 
bersome, since it does not rely on standard many-body 
techniques for fermionic or bosonic fields. 

An elegant derivation of the strong-coupling [Expan- 
sion has beeit-i pepently introduced by Sarkertj and 
Pairault et a/J23'Eil Metzner has obtained ±lie same re- 
sults, but his derivation is not as direct.E^I The basic 
idea of this approach is to decouple the intersite hopping 
term by means oLs. Grassmannian Hubbard-Stratonovich 
transformation.c^The resulting action for the Grassman- 
nian auxiliary field allows for a systematic perturbative 
expansion based on Wick's theorem. The (bare) propa- 
gator of the Tp field being proportional to the intersite 
hopping amplitude, a diagram with p lines (i.e. p propa- 
gators) is of order f . The Green's functions of the orig- 
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inal fermions are simply related to those of the tp field 
and can also be obtained in a systematic way. One of the 
difficulties of the strong-coupling perturbation theory is 
that there is now an infinite number of vertices. They 
are given by the connected atomic Green's functions and 
have therefore a non-trivial time dependence. High-order 
contributions are not easily obtained analytically, but can 
be nevertheless systematically computed using a specially 
designed computer program (see Ref. 20 ) . Another prob- 
lem lies in the fact that the calculated Green's functions 
(to a given order in t) do not have the correct analytical 
properties. This difficulty was circumvented in Ref. |2^ 
by representing the Green's function as a finite Jacobi 
continued fraction, which ensures the correct analytical 
behavior. 

Besides technical difficulties, the expansion around the 
atomic limit faces a more serious problem, namely the ex- 
istence of two dimensionless expansion parameters, t/U 
and t/T {T = 1/(3 is the temperature). Whereas hop- 
ping processes that do change the total number of doubly 
occupied sites are of order t/U, there are also processes 
that do not involve the small parameter t/U . The former 
processes correspond to virtual transitions between the 
lower (LHB) and upper (UHB) Hubbard bands, and the 
latter to intraband propagation. [Here we assume that 
the system is a Mott-Hubbard insulator.] This implies 
that expanding around the atomic limit is valid only at 
high temperature (T ^ t) , when intraband hopping pro- 
cesses are suppressed by the small parameter t/T, while 
we are mainly interested in the low-temperature regime 
T -^t. This failure of the strong-coupling approach was 
recognized early on, which motivated a reorganization of 
the perturbative expansion. Harris and Lange were the 
first to derive an effective Hamiltonian which describes 
exactly intraband hopping processes, while virtual inter- 
band transitianS|-ai'c considered perturbatively to a given 
order in l3 To first order in t/U, this yields the 

Hamiltonian of the t-J model. E3 

The aim of this paper is to modify the strong- 
coupling expansion for the Hubbard model introduced in 
Refs. [is] , ^ in order to clearly distinguish between intra 
and interband processes, thus allowing for an expansion 
in t/U. Hopping processes that do and those that do not 
conserve the total number of doubly occupied sites are 
identified by introducing slave particles (bosons) which 
act as projectors onto the empty, singly r^cupied, and 
doubly occupied atomic states (Sec. ||).E3 We empha- 
size however that the slave particles are treated exactly, 
which is possible since they intervene in fine only in the 
atomic limit. The Grassmannian Hubbard-Stratonovich 
transformation then requires the introduction of a two- 
component fermionic auxiliary field ip = , 
where ip+ (ip-) corresponds to a particle in the UHB 
(LHB). The single-particle propagator is naturally ex- 
pressed as a function of a 2 x 2 matrix self-energy. A 
simple approximation for the self-energy in the large-?/ 
limit predicts a metal-insulator transition. We compare 
this result with the one obtained within the Hubbard-I 



approxim ation . 



In Sec. HI 



we further modify our formalism in or- 
der to set up a strong-coupling expansion around a 
"non-degenerate" ground-state where each singly occu- 
pied state has a well defined spin direction (which may 
fiuctuate in time). This is achieved by introducing a 
time- and space- fiuctuati ng spi n-quantization axis in the 
functional integral (Sec. HI A). This provides a spin- 
rotation-invariant slave-boson description, which is an al- 
ternative-approach to the operator formalism derived by 
Li et al.tB The effective action of charge carr iers in the 
lower Hubbard band is derived in Sec. [II B. To low- 



est order in t/U, the coupling between holes and spin 
fiuctuations (which arise from the dynamics of the spin 
quantization axis) is described by a U(l) gauge-field the- 
ory, in agreement with previous conclusions To first 
order in t/U, we recover the action of thCj-^Xmodel in 
the spin-hole coherent-state path integraLEjO We also 
compare our results with those previously obtained by 
studying fluctuations around the large-[/ Hartree-Fock 
saddle point .Ej 



II. t/U EXPANSION 

We consider a £)-dimensional bipartite lattice. The 
Hubbard model is defined by the Hamiltonian 



H = -t ^ (c];^Cr'<T +h.C.) + JJ^nrT^ri, 
{r,r'),CT r 



(2.1) 



where c^a is a fermionic operator for a cr-spin particle at 
site r (cr =ti i)i ^ra = cl^Cra, and (r,r') denotes nearest 
neighbors. 

The hopping term in (2.1) does not distinguish be- 
tween processes that do or do not conserve the to- 
tal number of doubly occupied sites. Such a distinc- 
tion can be made by writing the fermion operator as 



(1 



f)Cr 



--1 (T) for a =T (i)]. 



It is more convenient to rewrite this decomposition using 
the slave boson representatioa of the Hubbard model due 
to Kotliar and Ruckenstein.Ea Introducing two fermionic 
(/r|, /rx) and four bosonic {er,Pr^ ,Pri, dr) operators, the 
four atomic states (empty, singly occupied, and doubly 
occupied) are expressed as 

|0,r) = eJlvac), 

1(7, r) =4,10, r) =C/rl|vac) ((T=t,i), 
I Ti, r) = ct^cijO, r) = dt/J^/Jjvac), (2.2) 

where |vac) denotes the vacuum of the enlarged Hilbert 
space. This enlarged space contains unphysical states 
which cari— |be eliminated by imposing the set of 
constraintsBa 

q(i) = ejer + Y,pLPra + dldv -1 = 0, 

Q^r} = fLha - pLPra - d^ = (a =?, i), (2.3) 



2 



at each lattice site r. The bosonic operators er, Pra 
and dr act respectively as projection operators onto the 
empty, singly occupied (with spin cr) a nd d oubly occu- 
pied states. The Hubbard Hamiltonian (2.1) is rewritten 
as 



where we use the notation 



H = -t J2 J2 ilLalr'c.'a+i^-C.) + Uj2dldr 
(r,r'),cr a,Q' — ± r 



where we have introduced the operators 



7r+<T = drplsfr. 



(2.4) 



(2.5) 



which annihilate a particle in the LHB an d th e UHB, 
respectively. The advantage of the form (2.4) of the 
Hamiltonian is that the hopping term now distinguishes 
between intraband (a = a') and interband (a = —a') 
processes. 

The partition function Z = Tre"'^'-^^^^-' can be writ- 
ten as a functional integral over fermionic (/) and bosonic 
{e,p,d) fields {N — J2ra fra frcr IS the_±otal number of 
particles and fi the chemical potential)La 

Z ^ J dX J 'D[f,e,p,d]exp^^-Sat[c,e,p,d; A] 

+ E / drrro.Jrr'lr'c.'a}, (2.6) 

^ r\j ' fr 



r,r' ,a,a' ,a 



where r is an imaginary time varying between and (3. 
trr' equals t if r and r' are first neighbors and vani shes 
otherwise. The variables "fraa are defined from (2^) by 
replacing the operators by fermionic or bosonic fields. 
A = (Ar^'jArCT^) denotes a set of time-independent La- 
grange multipliers that impose the constraints (|]^). Sa,t 
is the atomic part of the actionBj 

We are now in a position to apply the decoupling pro- 
cedure introduced in Refs. Sinc e the field jraa 
appearing in the hopping term [Eq. ( |2.6| )] carries a band 
index (a = ±), the auxiliary fermionic field involved in 
the Grassmannian Hubbard-Stratonovich transformation 
is a two-component field 



Ipr-a- 



(2.7) 



This is the main modification with respect to the strong- 
coupling expansion of Refs. The procedure to 
derive the action of the auxiliary field is similar and is 
briefly described below. The partition function reads 

Z = J V[^P]exp[-J2rJ-ti^b} 

a,b 

X dX / P[/,e,p,d]exp|-5at + E(^a7a + c.c.)|, 

J J a 

(2.8) 



'(pa = Ipr^aaCTaiTa): E= E j^Ta- (2.9) 



We denote by i^r' the hopping matrix to emphasize that 
it acts on the band index of the ip field. It can be written 
as the tensor product of irr' with a 2 x 2 matrix acting 
on the band indices: 



trr' 



1 1 
1 1 



(2.10) 



Note that the matrix t cannot be inverted. This is not 
a real difficulty since the final results can always be ex- 
pressed in terms of the matrix t, the inverse matrix 
appearing only at intermediate stages in the calculation. 
One could also consider a more general matrix where in- 
terband and intraband hopping processes have different 
amplitudes, which ensures the existence of the inverse 
matrix t~^. 

Performing the functional integration over the fields 
/, e,p, d and the Lagrange multipliers, we obtain 

fdX f P[/,e,p,d]exp{-5at + E('/'a7a + c.c.)} 

= Zatexp{w^[i^*,^]}, (2.11) 

where Zat is the partition function in the atomic limit 
(t = 0) and [■(/'*, ■(/'] the generating functional of the 
connected atomic Green's functions 



Thus the action of the tjj field is given by 



(2.12) 



(2.13) 



whereWhp* , Tp] can be obtained explicitly by inverting 
Eq. (|I|): 



= E E • • • V'LV'fcH • • • i'b.Gf:^ 



R=l 



(2.14) 



The primed summation in (2.14) reminds us that all the 
fields in a given product ■0ai ' ' ' "^fei share the same value 
of the site index. Note that the integration of the bosonic 
fields and the Lagrange multipliers has been done exactly. 
A difference with the strong-coupling expansion intro- 
duced in Refs. is that the effective action of the 



3 



auxiliary field ip involves the connected atomic Green's 
functions projected onto empty, singly occupied, or dou- 
bly occupied states. For instance, for the single-particle 
Green's function, we find (see Appendix ^ 



G+a{iuj) = 



1 



flrr 



(2.15) 



for the LHB and UHB atomic Green's functions, respec- 
tively. Here w is a fermionic Matsubara frequency, and 
ricr the average number of a-spin particles per site. The 
atomic Green's function of the original fermions (c) is 
simply given by the sum of G-a and G+o-- 

The free propagator of the auxiliary field is — t, i.e. 
— ('(/'oV'6)free = —tab, and wc may now use Wick's theo- 
rem to treat the interaction terms contained in ^^[-0*, ip] 
perturbatively. A diagram with p lines (i.e. p free prop- 
agators) is of order . The matrix structure of t allows 
to distinguish between processes of order t/U and those 
of order t/T. The former correspond to off-diagonal ele- 
ments of trr' and the lat ter t o diagonal elements. 

Although the action ( ^.13| ) allows to set up a system- 
atic t/U expansion, some difficulties remain. The first 
step would be to find an exact or approximate solution 
to lowest order in t/U , before considering interband cou- 
pling perturbatively. But this requires to consider an infi- 
nite number of vertices, with no obvious small expansion 
parameter since t/T ^ 1 at low temperature- Except in 
one dimension, where an exact solution existsjj this prob- 
lem is still unsolved. In Sec. Ill, we propose a slightly 



different formulation which solves (at least partially) this 
problem. Nevertheless, it is interesting to discuss approx- 
imate expressions of the single-particle Green's function 
that can be obtained from Eq. (2.13). This is done in the 
next section. 



Single-particle Green's function 

The single-particle Green's function Qab = ~{lall) is 
related to the propagator oL-the auxiliary field Tab = 
— {ip*ipb) by (in matrix form)L£l 



(2.16) 



Introducing the self-energy T of the auxiliary field, T ^ 
—t^^ — r, we obtain 



t + r- 



(2.17) 



The self-energy F is now a matrix with respect to the 
band indices. In reciproqual space, ro-(k, iw) is a 2 x 2 
matrix. Consider a i?-particle vertex for the field ^ 
[Eqs. ( ^.13D and ( |2.14| )]. U p {p < R) incoming parti- 
cles are in the LHB, then there are exactly p outgoing 
particles in the LHB. One can then easily convince one- 
self that the 2x2 matrix self-energy Fa (k, is diagonal 



(even though the hopping matrix t has off-diagonal ele- 
ments with respect to band indices). We denote its two 
components by r_cr(k, iw) and r+„{\i,iu}). 

The Green's function of the original fermions is related 
to g by 



a,a.' 



From Eq. (2.17), we obtain 

n°*(k,ic^) 



l-f-ikr^,°t(k,zt^)' 



where 



r^°*(k,zc^) = ^r„,(k,»c.). 



(2.18) 



(2.19) 



(2.20) 



and is the Fourier transform of irr' • [In the absence 
of interaction, the band energy is ek — ^ ^kj Eq. (2.2C) 
is nothing but the result obtained in Ref. |2y, where F'°* 
was calculated in perturbation theory with respect to the 
hopping amplitude t. To the extent where the self-energy 
F can be calculated exactly, our strong-coupling expan- 
sion is similar to that of Ref. PQ. 



The large-U limit 

Interband transitions are suppressed when J7 ^ oo, so 
that the off-diagonal elements of tk do not play any role. 
In the large-C/ limit, it is natural to make a diagonal ap- 
proximation where interband transitions are neglected. 
The Green's function is then determined by 



ik + F+o-(k, icj) 1 






tk + F_„(k, zcj) 



(2.21) 



Eq. (2.21) gives 



t/cr(k,ia;) 



+ 



1 



fk + F+,(k,zcj)-i ' ik + F_<,(k,ia;)-i' 

(2.22) 



The calculation of F+j and F_ct remains a very difficult 
task. The simplest approximation (lowest order in t) con- 
sists in retain ing o nly the Gaussian part of the action 



(2.23) 



From (2.22) we deduce 

g„(k,iu}) 



+ 



icj + ^ + (1 — ■ns)t]^ iuj + fj. — U + rigt)^ 

(2.24) 
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Retai ning only the Gaussian part of the action S[ip*,if:] 
[Eq. ( 2.23| )] is clearly unjustified, since this amounts to 
ignoring terms of order {t/Ty (p > 1) which are not 
small at low temperature. Nevertheless, this approxima- 
tion yields interesting results. At half-filling (/i — U/2 
and = 1/2), we obtain two Hubbard bands with dis- 
persions 



2 ' 
U + t^ 



(2.25) 



The Hubbard bands can be seen as due to the broad- 
ening of the atomic levels when the intersite hopping is 
switched on. They are separated by the Mott-Hubbard 
gap A = U — B /2 {B is the total bandwidth in the ab- 
sence of interaction), which closes for the critical value 
Uc = B/2 oi the Coulomb interaction. This simple 
approximation [Eq. ( p. 24 )] therefore predicts a metal- 
insulator transition at half-filling. It should be noted 
however that the bandwidth of the Hubbard bands is 
only half the total bandwidth B in the absence of inter- 
action. This is clearly an artifact of our approximation 
where only the particle (or the hole) that was added to 
the system can propagate. If the particle injected in the 
UHB has spin a, it can hop to the neighboring site only 
if the latter is occupied by a tr-spin particle, so that the 
hopping amplitude is reduced by a factor ng. [A similar 
results holds for the propagation of a hole in the LHB.] 
We expect the exact bandwidth of the Hubbard basds to 
be of order B in agreement with Mott's arguments. c3 For 
the same reason, we find that the Mott gap A separat- 
ing the two bands closes for the critical value Uc = B/2, 
while Mott suggested Uc ^ B. 

The fact that the single-particle propagator is natu- 
rally expressed as a function of a 2 x 2 matrix [Self-energy 
has also been noted by Logan and Nozieres.tJ This re- 
sult is reported in Ref. ^ where an expression similar to 
Eq. is discussed [Eq. (3.38) of Ref. |]. The differ- 

ence with Eq. ( 2.24 ) is that the Hubbard bands have the 
full bandwidth B. 

The result obtained in this section should be com- 
pared with the Hubbard-I approximation,^ which pre- 
dicts the opening of a charge gap for any finite value of 
the Coulomb repulsion U . The Hubbard-I approximation 
is obtained from Eq. (2.19) with the lowest-order contcLi 



bution in t for the self-energy r'°*, i.e. r^°* = y^^Ga^.El 



It can also b e obt ained from Eqs. ( p. 21 ) and (2.23) by 
including in ( 2.21 ) the off-diagonal elements of the ma- 
trix ik- Thus we see that the Hubbard-I approximation 
includes interband coupling without giving a correct de- 
scription of the Hubbard bands to zeroth order in t/U . 



III. NON-DEGENERATE t/U EXPANSION 

An important feature of the strong-coupling perturba- 
tive theory is that one expands around a state which has 



a huge degeneracy near half-filling, since a particle on a 
singly-occupied site can have its spin up or down. In this 
section we propose an alternative approach which allows 
to expand around a "non-degenerate" ground-state. 

In our formalism, the degeneracy of the ground-state 
shows up in the atomic Green's functions G^'^ which con- 
tain all the information about the atomic limit {t = 0). 
These Green's functions have to be calculated with the 
chemical potential /j, defined by the full Hubbard model. 
At half-filling, ji = U /2 due to particle-hole symmetry. In 
the presence of a small concentration of holes, the chem- 
ical potential lies near the top of the LHB, i.e. above the 
energy of the isolated atomic state. Therefore, at low 
enough temperature (when T is much smaller than half 
the bandwidth of the LHB), the system in the atomic 
limit is in its ground-state with exactly one particle per 
site (even in the presence of a finite concentration of 
holes). To suppress the degeneracy of the ground-state, 
we then proceed as follows. We impose the particles to 
have spin up. In order to restore spin-rotation invariance, 
we then allow the spin-quantization axis to fiuctuate in 
time and space, and integrate over all possible configura- 
tions in the functional integral. For a given configuration 
of the spin-quantization axis, the expansion in the hop- 
ping amplitude t becomes a non-degenerate perturbation 
theory. We shall show that this formulation allows to 
derive the effective action of charge carriers in the LHB 
to first order in t/U . We recover the action of the^-pZ. 
model in the spin-hole coherent-state path integralHsS 
We also comment on the effective action obtained by 
studying fiuctuations around the large-[/ Hartree-Fock 
saddle point Within our formalism, the results of this 
approach are reproduced by retaining only the quadratic 
part of the auxiliary-field action. This Gaussian approx- 
imation however does not capture all processes of order 
t/U , and therefore does not allow to derive the t-J model. 

In section HI A , we express the partition function as a 
functional integral where the spin-quantization axis fiuc- 
tuates in time and space. This provides a spin-rotation- 
invariant slave-boson approach to the Hubbard model. 
An alternative approach, based on-lhe operator formal- 
ism, has been proposed by Li et alxn The effective action 



of charge carriers in the LHB is then derived in Sec. IIIB 



A. Spin-rotation-invariant slave-boson approach 

The procedure to introduce a fluctuating spin- 
quantizati nn , ax is in the functional integral has been given 
by Schulz.Eacj It simply amounts to introducing a new 
field related to the old one by a unitary matrix i?r, which 
rotates the spin-quantization axis at site r and time r 
from z to a new axis defined by the unit vector Jlr (z 
is the unit vector along the z axis). In the present case, 
the procedure is slightly more complicated due to the 
presence of the slave bosons. Indeed it has been shown 
in Ref. that the bosonic field transforms under a 
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spin rotation as a 2 x 2 spin matrix field Praa' rather 
than as a two-component spinor field p^^. We show be- 
low that the introduction of a 2 x 2 spin matrix field can 
be avoided if one builds the functional integral from the 
very beginning. 

Let us first go back to the operator formalism. At each 
lattice site, we rotate the local spin-quantization axis by 
introducing a new fermionic operator 



where Cr = (Crt,Crx)^, (^r = (0rT,0ri 

unitary SU(2)/U(1) rotation matrix defined by 
RrCJzRl = cr • J7r. 



(3.1) 

)'^, and i?r is a 



(3.2) 



= (cajja^jCT^) stands for the Pauli matrices. The 
field (/>r has its spin-quantization axis along the (time- 
independent) unit vector Jlf. The U(l) gauge freedom 
is due to rotations around the z axis that do not change 
the state of the system. A convenient gauge choice is 



Rr 



cos(6'r/2) -e-'"^-- sin(6ir/2) 
e^^- sin(6'r/2) cos(6lr/2) 



(3.3) 



where Or,(pr are the usual polar angles determining the 
direction of fir- 

The interaction term being rotation invariant, the 
Hamiltonian is rewritten as 

H ^ - Y,4>lRUrr'Rr4r' +UY,4'l^4>r^4'll$rl. (3.4) 



As in Sec. we introduce fermionic (fra) and bosonic 
(er,Prcr,dr) Operators at each lattice site (with spin- 
quantization axis given by fir) which satisfy the con- 
straints (2.3). In this slave boson representation, the 
Hamiltonian is given by 



"time" Tk = kp/M (fc = 0, • • • M) and therefore becomes 
a dynamical variable ft^^T) in the continuum time limit 
(Af — > oo). Spin- rotation invariance is obtained by sum- 
ming over all possible configurations of the unit vector 
field rJr, i.e. 



r k=l 



vet z. 



(3.7) 



The details of this derivation are given in Appendix 
We obtain 



Z = J Vfl J dX J V[f,e,p,d]e- 



(3.8) 



where A = (Ar^^Ar^) denotes a set of time-independent 
Lagrange multipliers that impose the constraints ( p.3| ). 
The action is given by 

S = S,t- J2 j d^lLRltrr'Rr'lr'a'. (3.9) 
r,r',a,a' 

The atomic part reads 

'S'at = 5^?^ + J dr f*„p*^{RlRr)aa'fva'Pva', 



(3.10) 



where i?r = drRr- •S'at'' is the atomic actio n for a time- 



independent spin-quantization axis [see Eq. (B19)]. Fluc- 
tuations of the spin-quantization axis modify the hopping 
term [Eq. ( |3.9| )] an d ind uce an additional term in the 
atomic action [Eq. ( 3.10D ]. As shown below, the latter is 
related to the Berry phase term of a singly occupied site 
(which behaves as a spin). 



H 



E 

r.r' .a.a' 



Rv'lv 



U J2 dldr 



(3.5) 



B. Effective action of charge carriers in the LHB 



where = (>«!' 7rai)^- Iraa is defined by Eq. (|2J). 

In order to express the partition function as a func- 
tional integral, we divide the "time" interval (3 into M 
steps: 



Z = Tr( e"'(^"^^) • • • g-^^^-''^) 



(3.6) 



where e = (3/M . We then introduce M — 1 times 
the closure relation using mixed fermion-boson coherent 
states. An important point here is that the local spin- 
quantization axis depends a priori on the discrete 



J 



S — S'at ^ 



j dr-ilRlUr'Rv 



1 



The partition function ( |3.8| ) obviously preserves spin- 
rotation invariance. As discussed above, we can now as- 
sume that the LHB is populated only with up-spin p. 
tides (in the local spin reference frame defined by fir) 
This implies that only down-spin particles can be intro- 
duced in the UHB. Spin- rotation invariance is maintained 
since we integrate over all configurations of the unit vec- 
tor field in the functional integral. 

The assumption that the LHB is populated only with 
up-spin particles is imposed by eliminating the pi field. 
The action then reads 
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•' r 

r cr 



(3.11) 



where 



7rT 



7r-T 
7r+i 



(3.12) 



The fields 7r-x ^-nd 7r+T do not appear any more. Note 
that the UHB (LHB) can be labeled by the band index 
a = + (a = — ) or the spin index cr =| (cr —^)- In 
the following we use the spin index to label the Hubbard 
bands. A^^ — (i?Ji?r)|f is defined by writing the gauge 
field AO = i?t4 as AO = 

The Hubbard-Stratonovich decoupling of the hopping 
term requires the introduction of the auxiliary fermionic 
field 



The partition function can be written as 

Z = j Vndet{i) J V[i^]e^p{~Y.^aiab^b} 

a,b 

X dX I?[/,e,p,d]exp{-5at+E(V^7a + c.c.)}, 

(3.14) 



where the hopping matrix 



(3.15) 



(3.13) 



acts on the spin/band indices. Note that det(£) depends 
on Jl and should therefore be kept explicitly in the func- 
tional integral. Using Eq. (3_^), the matrix R\Rr' can be 
expressed as 



J 



^t^_^/ (1 + Or-OrO^/V*^"/' (l-nr-nr')^/V<^'/2 
V2 I -(l-r2r-r2r')^^ 



'^^y'^ (l+nr-f^r')'/'e-**-'/2 



(3.16) 



r 



where 'I'r.r' = $(r2r,r2r') is the signed solid an- 
gle spanned by the vectors r2r,J^r',z, and ^'^ = 
<i>(r2r, — f2r') ~ 2(/3r'. $ and <i>' satisfy the relations 



+ 27r. 



Performing the integration over the fields f,e,p,d and 
the Lagrange multipliers, we obtain 



^ J at 

= Y.fdTAl^. (3.19) 



z = prjZat[rj]det(£) 



Higher-order corrections, of order t/U ^ are ignored. This 
is justified since the spin dynamics is slow compared 
to the interband-transition dynamics. [The second line 
, , of Eq. ( 3.19| ) is easily obtained by following the "phi- 

I?[i/']exp|-EV'aC6Vfc + W^[V'*,V';f^]}, (3.17) losophy" of Appendix 1^.] 5b is a collection of ^Berry 



a,h 



where 



Zat[f2]= j d\ j V[f,e,p,d]e^v[~S^t] 

){-5b[0]} (3.18) 



phase terms for spins localized at the lattice sites.l 
ing Eq. (B.3), it can be written as 



SB[n]^'-J2jdT{l- cos 0r)iPr 



Us- 



.(0), 



E J drAifl) ■ n, 



(3.20) 



is the partition function in the atomic limit for a given 
configuration of the fluctuating spin-quantization axis 
Clr- The generating functional of the connected atomic 
Green's functions Wlip* depends on Cl^ since S'at 
does. ^'^[Jl] is calculate d by treating the gauge field 
in perturbation in Eq. ( |3.18| ). To lowest order (in a cu- 
mulant expansion) 



where A(rj) is the vector potential created by a magnetic 
monopole sitting at the center of a unit sphere. 

Th us, th e acti on of the ip field can be written as [see 
Eqs. ( |3T7| ) and ( p^ j 

S'[V'*, V; f^] = - Indet(t) + 5b [fl] 



7 



a.b 

- lndet(£) + Sein] + ^ raiU + Gab) A 

a,b 

i E Gi\a„f„5.cv':>f>>b,. (3.21) 

ai,bi 



In the second line of (3.21), we have retained only the 
quadratic and quartic parts, since this is sufhcient to ob- 
tain the effective action of holes in the LHB to first order 
in t/U (see below). 

In order to completely determine the action 
S'[^/'*, ■0; we need to compute G and G^^"^. The single- 
particle Green's function G is calculated in Appendix 



0: 



= - sgn(a)^! 





rz ' 



(3.22) 



where sgn(f7) = 1(— 1) for a =1 (|). G'^"-' is the atomic 



Green's function corresponding to 5*^^ [Eq. ( |3.11 )] 



Gf\T)=ei-r + rj)e^\ 
Gj°^(r) = -0(r-77)e(^~^)", 



(3.23) 



where 77 ^ 0+, and the limit T ^ has been taken. In 



Fourier space, Eqs. ( 3.23 ) become G'"^\iuj) = {ito + ^j.) 



and Gf\iuj) ^ {ioj + - U)-\ 

Neglecting the Berry phase term A^, we can obtain 
G^^° using the method of Appendix It is straight- 
forward to modify this approach in order to impose the 
constraint that the LHB is populated by up-spin particles 
only. At T = 0, we find 



||(Ti,r2;T3,T4) 



= 0, (3.24) 

= G\°\t, - Ts)Gf\T2 - n) 

x[9{n - T2) - 9{t3 ^ Ti)]. (3.25) 



Since Gj/^ — 0, G^^^ necessarily involves (virtual) in- 
terband transitions. The typical energy scale (~ U) for 
these transitions being much larger than the typical en- 
ergy scale for spin fluctuations (^ J = 4t'^/U), we can 
Ignore the effect of A° on G"^ Eq. ( |3.24| ) can easily be 
extended to higher-order connected Green's functions: 



G 



Rc 



(i? > 2) 



(3.26) 



In the limit U 00, where interband transitions are 
suppressed, the action of the ip field is Gaussian. It re- 
mains nevertheless highly non-trivial since the intersite 
hopping term depends on the spin variables fir- [This 
limit is discussed in Sec. [IIBl.] Furthermore, this im- 



plies that only a limited number of higher-order vertices 
G^'^ {R > 2) is necessary to obtain the effective action to 
a given order in t/U (if one considers only the quadratic 



and quartic parts of the effective action) , since these ver- 
tices necessarily involve interband transitions. In partic- 
ular, it is sufficient to consider G^^"^ to obtain the effective 
action to first order in t/U. This is an important differ- 
ence with the approach of Sec. where the action of the 
auxiliary fermionic field ijj involves an infinite number of 
vertices even in th e limit L/ — > cx). 

From Eq. ( 3.21 ), it is now possible to integrate out 
the UHB to obtain the effective action of a hole in 
the LHB. It is more convenient to first perform an- 
other Hubbar d-Str atonovich transformation of the hop- 
ping term in (3.21). If this transformation were carried 
out exactly (i.e. on the exact action S[^p* ,tp;0,]), one 
would essentially "undo" the first Hubbard-Stratonovich 
transformation and recover the original action S['y* , 7; fl] 
[Eq. (3.11)]. Here the idea is that to first order i n t/U , 
one can truncate the action S[ip*,tp; fl] as in Eq. ( 3.21 ). 
One then obtains an effective action S[^*,j;n] which 
contains all processes of order t/U. It is then possible to 
integrate out the UHB (i.e. the field 74) to obtain the 
effective action S'[7p7|;n] for the LHB. 

Following this procedure, we obtain 

JvnJ I?[7] exp{-5[7*,7;fl]}, 
5'[7*,7;^^] = Sb[^] -^lliablb 

a,b 

-W^[7*,7;ll] -lnZ[n], (3.27) 



where 



z[n] 



V 



[ib]e,^p{w[r,^;^]}- (3.28) 



Note that 7, which is the auxiliary Grassmannian field 
of the Hubbard-Stratonovich transformation, is now an 
independent variable. Vl^[7*, 7; fl] is the generating func- 
tional of connected Green's functions obtained from the 
action -W[tp* ,tp;fl\. To lowest order in t/U, lnZ[n] = 
Indet(G) = S'B[f2]. Again we neglect corrections of or- 
der t/U to the Berry phase term. Retaining only the 
quadratic and quartic parts of the effective action, we 
have 



S[r,TM = -Y.^:{iab^Gab)lb 



a.b 



"(2ff E ^"a2,&i&2 7ai7a2 7b2 7fci, (3.29) 



where 



Gab 



'{Ml) 



-w, 



(3.30) 



In Eq. (3.29), it is sufficient to determine G and G^^"^ to 
first order in t/U . We find 



ab — ab + ^ 



ab ' 



G 



lie 

0102,6162 



'-'aia'j'-'a24'^a;4,6;6^^6'j6i'-'6^62' y^'"^) 



a'b'. 
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where 



(3.32) 



a' ,6'. 



is the "one-loop" correction to —G ^. From (3.31), we 
conclude that G^^'^ is simply related to the atomic two- 
particle vertex: 



(3.33) 



We thus obtain 



+ ]W ^ r;,\a2.fti62 7ai7a2 7&2 76i- (3-34) 

Note the important role played by If we consider 

the one-loop correction due to T^^, we obtain (among 
other contributions) an atomic contribution to the prop- 
agator of the field 7, which is already included in G. 
This contri bution is precisely canceled by (see sec- 
tion 

By per- 



[IIB2I) 

Eq. (3.34) has a clear physical meaning. 



forming two successive Hubbard-Stratonovich transfor- 
mations, and truncating the action S[ip* , ip; 17] in the in- 
termediate step, we have effectively summed the atomic 
contributions to the single-particle propagator and the 
two-particle vertex. Since Fj/^ — 0, perturbative cor- 
rections (due to Fj^j^ ij^) to the Gaussian action are of 
order t/U. In the limit [/ ^ 00, S'[7*, 7; fJ] reduces to its 
Gaussian part (except the term 7*G'(^)7 which is of order 
t/U). The effective action for a hole in the LHB can be 
derived by integrating out the UHB, the two-particle ver- 
tex F|| II being considered within perturbation theory. 



1. Limit U ^ 00 



where we have used Eqs. (3.16), (3.22) and (3.23). From 
now on, we drop the spin index for fermions in the LHB 
and denote by 7 the field 71 . The action ( 3.35| ) is similar 
to that obtained by, Schulz from the large-J7 Hartree- 
Fock saddle point .£3 At half-filling (7r|7rT = l)j it re- 
duces to 5*5 [rj]. As expected, the half-filled Hubbard 
model becomes a collection of independent spins in the 
absence of interband coupling (U — + 00). The coupling 
between a hole in the LHB and spiiL. fluctuations is de- 
scribed by a U(l) gauge- field theory.E3 t3 As pointed out 
in Ref. the introduction of a hole (absence of a 7| 
particle) in the LHB has two effects, (i) The Berry phase 
term A'^^^*^^, which is alive when the site r is occu- 
pied by a particle, is suppressed by the introduction of 
a hole. This conclusion was first obtaiped by Shankar 
from semi-phenomenological arguments Hil (ii) The hole 
can propagate through the system, but spin fluctuations 
modify the kinetic energy term. The hopping term be- 
tween two neighboring sites r and r' has its amplitude 
reduced by a factor [(1 + fl^ ■ flr')/2]^^'^ , and acquires 
the phase <i>r.r'/2- The latter can be interpreted as re- 
sulting from an effective magnetic field depending on the 
spin configuration.^ It is then equal to the circulation of 
the magnetic vector potential on the link between r and 
r'. When going around an elementary plaquette, the ac- 
cumulated phase <&i234 corresponds to the magnetic flux 
through the plaquette, and is given by half the solid angle 
spanned by the unit vectors Jli, 172, J^s, 174: 



1 



1234 



[*1 



2,3 



3,4 



<i>4 



(3.36) 



We expect the kinetic energy term to be optimized for 
a ferromagnetic configuration of the spins, since in that 
case the effective magnetic field vanishes and the hopping 
amplitude takes on its maximum value|_.Xlpis, is nothing 
but the familiar Nagaoka phenomenon O'HH 



In the absence of interband coupling, we obtain the 
following effective action for the LHB 

^LHB=E / rfr7;(5.-A^ + A°j7r 

r 

'7f E / rfr(l + 17r.l7rO^/^[7;7r'e^*-'/' + c.c.], 

(r,r'> 



2. Correction of order t/U 



The integration of the field 7^ in Eq. (3.34) will gener- 



ate a correction S'l^b to the effective action S^^^ . Con- 



sider first the quadratic part of S'lhb- To first order in 

,(1q) Mlb) 
LHB ~r ^LHE' 



(3.35) t/U, it is given by S'l^b + '^'lhb' where 



^[hb= E /drdr'7;(r)(frr'(r))TiGi(r,r')(W'(r'))iT7r"(r'), 

r,r' ,r" 

'5'lHB E/ 0'TrfT'dT2dT4 7;(T)7r(r')F|\ ||(T,r2;T',T4)(G^^ + ^"u)rr4,rr2- 



(3.37) 
(3.38) 
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S'lhb '^^^ expanded in a series in powers of in- The 
lowest-order term exactly cancels the term 7*G'j^^7| in 
Eq. ( |3.34 ). The next-order term (0(i|j^)) is proportional 
to / rfT2rfT4 G||'^||(ti, r2; Ta, r4)Gj^ (r4, and therefore 
vanishes [since G|^j^^|j^(ri, r2; T3, T4) cx Gx(r2,T4)]. For 
a similar reason, all higher-order terms vanish. Thus 
the quadratic part of the effective action reduces to 



i.e. 7r (''')7r' ('''') — 7r(''')7r' (''')• The simplc replacement 
7r/(T') — > 7r'(''") would howcvcr lead to a wrong result. 

Instead, we write 7(t') — 7(r)e~'^''^~'^ which follows 
from the equation of motion for the field 7 in the atomic 
limit .El S'L^g is then approximated as 



5, 



(0) 
LHB 



(la) 
LHB- 



c,(la) _ 
•^LHB — 



Consider now S'^^g (shown diagrammatically in 
Fig. |l|(a)). We have argued above that we can ignore the 
effect of spin fluctuations on (virtual) interband transi- 
tions. This allows us to replace G|(r, r') by G^^\t — r') 



dr 7; (T)7r" (r)(<rr' {t))u iU'r" (t))xt 

J dr' Gf^ (t - T')e-^("-"') . (3.39) 



and {tr'r"{T'))u by (tr'r"(r))tx in Eq. ([3.371) . Further- 
more, since the characteristic energy scale (~ U) of a 
virtual transition is larger than any other typical energy 
in the system, we can assume S^^^ to be local in time. 



Using Eq. ( 3.23 ), the sum over r' can be written as 

- TrfT'e-^^"-"') ~ (3.40) 
Jo U 



Using Eq. ( 3.16 ), we finally deduce 



^LHB 



2U 



E / '^^ W7r + 7;'7r'](l - f^r • ^r') 
^ E /dr(l-J^r-f^rO'/'(l-f^r'-Or,0^/'[7;7r''e^'*-'-*^''-'^/'+C.C.], (3.41) 

Ir. r-l r.ll\ 



Here 
(r,r') and 



r,r',r") 



(r,r',r") ' 

stands for a three-site term with both 
' first neighbors (r 7^ r"). 



r 



The effective action 5, 



(0) 
LHB 



"^LKB i^ similar to that of 



Ref. |35| derived by studying fluctuations around the large- 
IJ Hartree-Fock saddle point. In the absence of holes in 
the LHB, it reduces to the action the AF Heisenberg 
model with coupling constant J = 4i^/[/, 

5Hcis[0] = + J E /^^ ( "^ 4"'" - -X (3-42) 

which is the expected result. Nevertheless, we will show 



(0) 
LHB 



■^LiJB ^''^ effective action to first or- 



that 

der in tjV . One has to consider the quartic contribu- 



tion iSlhb to S'lhb which results from the integration of 



^(1) 



the UHB. It will turn out that "S'lhb' when written in a 
time-ordered fashion (as it should be in the functional in- 
tegral), generates a quadratic term which exactly cancels 
"^LHB- What will be left is nothing but the action of the 
t-J model. In fact, it is easily seen on physical grounds 
that iS'L^g is not the action that we expect to first order 
in t/[/. The introduction of a hole at a given site r should 
suppress the AF exchange with all the neighboring sites. 
Now, according to S'lhB' ^^e AF coupling between two 
neighboring sites r and r' is reduced by the presence of 
a hole at site r or r', but vanishes only if two holes, one 
at each site, are introduced in the system. 

The quartic contribution to the action, shown in 
Fig. |^(b), is given by 



^mB= E / dndT^dr^dn J dT^dTiT]l^^{n,r^;T3,T:,)Gf\r2~T^)Gf\T:,-n) 

r.r'.r" 

x(fr'r(T2))Ti(trr"(T4))iT7^(n)7^'(T2)7r"(T4)7r(T3), 



where 



r^\,T^(ri,r2;T3,T4) = G;°)-'(ti - r3)Gf^-'(r2 - nMn - r^) - 9in ~ t^)\ 



(o)-i. 



(3.43) 



(3.44) 



As for S'l^3, one can assume S'lhB be local in time. Using 7(r) ~ 7(0)e^'^ and 7*(t) ~ 7*(0)e we obtain 

•^[1115= E / ^^l7r(n)7r'(n)7r"(Tl)7r(n +?7)(tr'r(Tl))Ti(<rr"(Tl))xT 
r,r' ,r" 

X J dT2dT3dne-''^^-+^--^'-^'^Gf^~\n - r3)Gf{T2 ~ t^)\B{t^ - r^) ~ 0{t3 - t^)], (3.45) 
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where ry — s- 0"*". This equation has been obtained by noting that 

Gf ^'^T) = fiSir + 7]) - 5{t + Tj), (3.46) 

where the overdot denotes time derivative. The infinitesimal 77, which is required by a proper treatment of equal-time 
correlation functions, is crucial here. It implies that T3 is larger than ti in Eq. (3.43). When approximating S^^^ by 
a local vertex (in time space), one should keep track of this time ordering, since it do es no t correspond to the one 
implicitly assumed in the functional integral. Performing the sum over T2,T3,r4 in Eq. (3.45), we deduce 

^hI--^ E /^^7;(r)7r(r + r;)7rMr)7r"(T)(4'r(r))T,(4r"(r))iT. (3-47) 

r,r' .r" 

Noting that Jri'^)lr{T + rj) = 7r (''')7r(''') — 1, we eventually come to 



g(lc) _ _ Mia) 
'-'LHB — "^LHB jj 



^ E / dri;iri;nr"{U'r)niirr")lV (3.48) 



The quadratic contribution which is generated by the correct time-ordering in the functional integral cancels S^^^^. 
The effective action to first order in t/U, S'lhb + "^lhb + '^'lhB' therefore given by 



5lHB = E / ^^7r (5. - A* + ^rj7r - ^ E l dr [l + ■ n^^ f'^Yll^r'C^"^-' 



/2 



+2 



^- E /dr(l-flr-J^rO'/'(l-f^r-Or.O^/'[e^^<^'-*"''^/^;7r7;'7r''+C.a^ (3.49) 



2[/ 

{y'y.v 



The action (|3.49|) exhibits the correct physical properties. T he Berr y phase term and the modification of the kinetic 
energy term by spin fluctuations have been discussed in Sec. IIIB 1. The J-term describes AF exchange interactions 
between neighboring sites (r,r'). It vanishes as soon as a hole is present at site r or r'. The three-site term is also 
familiar from the derivation of the t-J model. 



3. Relation with other formalisms 



In this section, we show that-the action ( 3.49| ) corresponds to the action of the t-J model obtained in the spin- hole 
coherent-state path integral.tSEZi We then make the connection with the slave-fer mion f ormalism. 
We first perform a particle-hole transformation: 7r ^ /ij, 7* ^ hr- The action (3.49) then becomes 



^LHB = Y. I dr [K[dr - tlh)K + Al,{l - KK)] 
r •' 

+ 7f E J dr {1 + fir ■ ^r'Y^^lKh, 
{r,r') 

+^E fdr[^^^^\){i-KK){i-h:,K,) 



e-*.,.'/2 + c.c.- 



(r',r,r") 

x[e'(*:..'-*:..")/2(i _ hlK)hl„K, + c.c.]. 



(3.50) 
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where /i^ = ~ /i is the chemical potential of the holes. 



The action ( 3.5C ) is precisely the action of theji-iJ-jnodel 
in the spin-hole coherent-state path integral.E30 Note 
that the three-site term (the so-cal led "pair-hopping" 
term, last term of the rhs of Eq. ( 3.50D ) is usually omitted 
in the t-J model. 

The rotation matrices Rr are elements of SU(2)/U(1). 
Instead of integrating over this manifold in the func- 
tional integral, one can choose to integrate over all 
SU(2) matriceg-|(with the properly normalized integra- 
tion measure) ,E3 i.e. one writes 



1 



(3.53) 



Note that we have written A'^ in a symmetric way. In- 
troducing new variables (Schwinger bosons) defined by 



Zr. 1 



-hlh. 



Rr - 



Zir 
Z2r 



with the constraints 

\4r\ + \4r\ - 1- 

Writing = (zir, •Z2r)^, one then has 



(3.51) 
(3.52) 



the constraints ( |3.52 ) become 
and the action reads 



Khr = 1, 



(3.54) 



(3.55) 



J 



^LHB = E / ^l^ridr - fih)hr + ^(fo^&r ^ b%)] + t ^ fdr [blKKK, + C.C.] 



E / dr [blabr ■ blabr' - (1 - h;hr){l - K,hr')], 

(r,r'> 



{3M) 



where, for simplicity, we 
hopping" term. The action 



have omitted the "pair- 
(3.56), together with the 



constraints ( 3.55| ), is the action of the t-J model in the 
slave-fermion formalism [see, for instance, Ref. 



IV. CONCLUSION 



We have reconsidered the strong-coupling expansion 
for the Hubbard model introduced by Sarkerllj and 
Pairault et alEB The main modification is twofold, (i) We 
introduce two different fermionic fields, 7+ and 7_, which 
correspond to particles in the UHB and LHB, respec- 
tively (Sec. H). (u) We organize the strong-coupling ex- 
pansion around a "non-degenerate" ground-state where 
each singly-occupied atomic state has a well defined spin 
direction (which may fiuctuate in time) (Sec. III). The 
action of the system, S'[7*,7;r2], is then expressed in 
terms of two fermionic fields, 7_| and 7+1 , and a unit vec- 
tor field Jlr- 7-T &n.d 7+1 correspond to particles in the 
LHB and UHB, respectively, the local spin-quantization 
axis being specified by the time-fluctuating unit vec- 
tor field fir- As in Ref. ^ the strong-coupling ex- 
pansion is carried out by introducing a fermionic auxil- 
iary field ip = ("0-11 V'-i-i)"^ s-nd performing a Grassman- 
nian Hubbard-Stratonovich transformation of the inter- 
site hopping term. 

Compared to the formalism discussed in Ref. our 
work brings two major improvements, (i) The expan- 
sion involves a single dimensionless parameter, t/U, and 



therefore does not break down at low temperature. In 
the limit U — > 00, the action S'[7*,7;r2] becomes Gaus- 
At finite 



sian (Sec. [II Bl). At finite the corrections to the 
Gaussian action can be obtained perturbatively by con- 
sidering a finite number of vertices, (ii) By introduc- 
ing a fluctuating spin-quantization axis, we have effec- 
tively suppressed the huge spin degeneracy of the atomic 
ground-state. Indeed, for a given conflguration of fir, 
we now expand around a "non-degenerate" atomic state, 
where each singly-occupied site carries a well defined spin 

(f2r/2). 



This formalism allowed us to recover the action of tJiei^ 
J model in the spin- hole coherent-state path integralHj'EZl 
It can also be used to directly study the Hubbard model. 
It is clear that the main technical difficulty comes from 
the spin variables fir, since the fermionic fields 7_| and 
7+1 can be integrated out in a systematic t/U expan- 
sion. Thus further progress is tied to an approximate 
treatment of spin fluctuations. The simplest approach 
consists in expanding around a broken-symmetry ground- 
state by making a saddle-point approximation on the spin 
variables fir- In the half-filled Hubbard model, a natu- 
ral choice is fir ~ (— l)'"z, which corresponds to an AF 
ground-state. Other choices, such as ferromagnetic or 
spiral orders, are also possible. Work along these lines 
will be reported elsewhereJlJ 
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APPENDIX A: 



In this Appendix, we show how we can obtain the atomic Green's functions projected on empty, singly occupied, 
and doubly occupied sites. One could directly use the slave boson formulation by explicitly integrating over the fields 
f,e,p,d and the Lagrange multipliers A. This is the method used in Appendix ^ Here we present a different (and 
more direct) approach. We consider a single site and drop the site index. 

In the basis {|0), | t), | i), | Ti)}: the operators jaa are given by (in matrix form): 



7+T 



7-T 



From the equations of motion, we then deduce 
where ji^}{T) ^ U {-T)^'i^lu [t) , and 
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(Al) 



7l,(r)=e-''^7L, 



U{t) = e 



/ 1 

e''^ 

e^^ 

y e^^A'-t^)^ 



(A2) 



(A3) 



Using (Al) and (A3), it is straightforward to calculate the Green's functions. For the single-particle Green's 
function, we find (for r > 0) 



G-.(r) - --TT[Um-.{r)llAQ)] 
z 

G+Ar) = --TT[Um+ArW+AO)] 



where 

z = Tr f/(/3) = 1 + 2e'^'' + e 
is the partition function of a single site. In Fourier space, we obtain 

G_<j(lLj) = 



P{2ti-U) 



(A4) 
(A5) 



G+c(iLj) = 



riif 



where 



iuj + j.1 — U' 



(A6) 



(A7) 



is the mean number of cr-spin particles per site. The method is straightforwardly extended to the calculation of 
higher-order Green's functions. 
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APPENDIX B: 



In this Appendix we derive the spin-rotation-invariant action given by Eqs. (3.9) and (3.10). 

In order to derive the functional integral, we introduce the mixed fermion-boson coherent states |C) = \f,e,p,d) 
defined by 



10 = expj^ [- frafL + ere]: + "^PrapL + dJl } | 



(Bl) 



where f^a is a Grassmann variable, and e^, Pra and are c-numbers. In Eq. (Bl), the spin-quantization axis is 
determined by the unit vector fir- The coherent states |C) satisfy the closure relation (/ is the unit operator) 



where 



drrfce-i«i'ic)(c|p-/, 



(B2) 



\Cf =Y[Y frafriy + eJCr ^ p^^Pra + d^rfr}, 
r a a 

lderdd*ddr Y\_{dfrcrdfradp*^dprcr)] , 



and A/" is a normalization constant. 



^=n(V'.onv;.o) 



(B3) 



(B4) 



is a projection operator which ensures that the closure relation acts only in the physical Hilbert space. The partition 
function is written as 



Z= / dC*rfCe"l'^l'(Cl^e-'^(^-''^)|0, 



(B5) 



where {(\ = (— /, e,p^. Splitting into M intervals of width e = P/M as in Eq. (3^), and introducing (M-1) times 
the closure relation (B2), we obtain 



/ lldadae-^^^^^'^^^'i[{a\p,e-<"-'^^^\a^, 

•' k=l k=l 



(B6) 



with the boundary conditions fra^M — —fra,o^ ^r,M — ^rS)i Pra,M — Prcr,Oi and (ir,M = c?r,o- The spin-quantization 
axis at a given site, defined by Rr,k<^zRl fc = ' ^r,k, may depend on the "time" k (but satisfies the boundary 
condition fir, a/ = J^r.o)- The projection operator P^. can be written as 



2tt 



r.k -*eA(l'Q<\' 

: p T,k^r,k 



n 



2tt 



(B7) 



Note that the operators Pk, Q^^l, and Q'^J are defined with respect to the spin-quantization axis flr^k- Since Q^^}. 

(2) . ' . . ' 

and Qj.J j, commute with the Hamiltonian H — fiN, we have 

/ \{dCkdCkdXke-^^-^\^''\"\{{Ck\e-'''''\Qk-i), 

•' k=l k=l 

= / n^CfcdadA,.e-S^^il'^'=l'n{(Cfcia-i)e-^^''(^-'^''-)} 

J 1 7 1 



(B8) 



fc=i 



k=l 



in the limit e 0. We use the notations 
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/--/n(^n 



2tt 

(i)n(i) 



r (T 



(410 



k-i/ 



(B9) 



The only difference with the standard derivation of the functional integralcEl for a mixed fermion-boson system 
comes from the evaluation of the scalar product (CfclCfc-i)- If -Rr.fc = Rr,k-i at each lattice site, we have the usual 
result 



{Ck\Ck-l) = exp|E[E/;„fe/r<T,fc-l + e;f.er,k-l + 'YPrcr,kPra,k-l + d^^fcdr,fc-l }, 

r a 



(BIO) 



since the spin-quantization is the same at "times" k and k—l. In the more general situation where the spin-quantization 
may fluctuate in time, we rewrite the coherent state |Cfe) as 



ICfe) = n{'°'"-''=^r,fe -Yfr<T,kfL^kPr<y,kpl^^k " /rT,fc/ri,fe/r\,fc/ri,fcO'r,fc4,fc } l^ac) 
r cr 

= II{'^"-.A;|0,r) -Yfra,kPra,kW,r,k) - /rT,fc/ri,ferfr.fc | Ti,l")}: 



(Bll) 



using the constraints (2.3). |0,r), jcr, r, /c), and | ti,r) denote the empty, singly occupied, and doubly occupied states, 
respectively [sec Eq. ( |2.2D ]. The singly oc cupie d state |cr, r, A:) depends on the "time" k since the particle has its spin 
quantized along the ^-axis. From Eq. (Bll), we deduce 



{Ck\Ck~l) = Y[{e;f,er,k~^l+Yfra,kPra,kfra',k-lPra\k-l{(J,r,k\a',r,k- 1) 



a. a' 



+ frl,kfrT,kfr'\,k-lfrl,k~ld*^kdr,k-l}- 



(B12) 



When Or. A; — ilr,fc -i, (g", r, k\a' ,r, k — I) = S^a'- Eq. ( B12 ) is then equi valen t to Eq. ( BIO ). [This equivalence is due 
to the constraints (2.3).] For fl^.k — ^r.k-i, this allows to rewrite Eq. ([B12[) as 



(CfclCfc-l) = exp|E \^ fc/ra,fe-l + e;_fcer,fc-l +Yp*„^kPra,k^l + d;,fcdr,fc-l 
r cT a 

+ Yfra,kPra,kfra',k~lPra',k-l{{Rl^kRv,k-l)aa' - S^^a') }, (B13) 
o-.cr' 

where we have used (u, r, k\a' , r, k ~ 1) — {rI i^Rr,k-i)aa' ■ Thus we may write the partition function Z as 

^ = Yl dQdCkdXk exp|- ^ ^ \^ /r<T,fc(/r^,fc - fra,k-l) + ej,fc(er,fc - er,fc-l) 
k=l k=l r a 

+ YPraAPra,k - Pra,k-l) +d* f.{d^^k ^ dr,fc-l) 

M 

frcr,kPr<j,kfra' ,k-lPva' M-l{iRl^kRr,k~l)aa' " ^^t^ct') - £ ^ -ftTfc (Cfc , Cfe- 1 ) } • 



(B14) 



fc=l 



Spin- rotation invariance is obtained by summing over all possible configurations of the unit vector r2r,fc, i.e. 



Z 



M 



Yin^^z. 



r fc=l 



47r 



(B15) 



When taking the continuum Hmit, one has to evaluate Kk{Q,Ck) (instead of Kk{CkXk-i))- Since all the operators 
involved in that quantity are defined with respect to the same spin-quantization axis (one can choose the axis Jlr,fc 
to express the Hamiltonian H — fiN), Kk{Q,Ck) is readily calculated. We finally obtain 
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Z = JvnJvxJ V[f,e,p,d]. 



where the action S is given by 
S 



r,r',Q:,a' 



5". 



(0) 



dT-fl^RlUr'Rr'lr'a', 

r,a,ar' 

E / drpA«+E/;.(5r 

r •' a 



ra'Pra' j 



-^J2X'^J+U)dr 



(B16) 

(B17) 
(B18) 

(B19) 



The variables 7 are defined in Sec. [I]. It is well known that in the standard slave boson approach, the Lagrange 



The 



multipliers can be chosen to be time independent since the constraints are preserved under time evolution! 
action S [Eqs. (B17-B19)] differs from the standard action by the additional term proportional to R^R which comes 
from time fluctuations of the spin-quantization axis. This term conserves the total number of p bosons and therefore 
the constraints Qi^\ It is also clear that it lets the combination f*^frcr — PraPro- invariant, so that the constraints 
QrCT^ are also conserved. Consequently, we can replace the functional integral over A(r) in Eq. (^ ) by an integral 



over a set A = (Ar^'', Ar^ ) of time-independent Lagrange multipliers: 



VX 



dX : 



n 



- f3dX 



(1) 



27r 



n 



2tt 



(B20) 



Note that the introduction of time-dependent Lagrange multipliers at an intermediate stage is what allows a simple 
evaluation of Kk{Cki Cfc-i)- To avoid this complication, one would have to consider the p operators as 2 x 2 matrices 
in spin space. 



APPENDIX C: 

In this Appendix, we calculate the single-particle Green's functions c'^a^ (r) and Go- (r) corresponding to the action 
and 5'at, respectively [Eqs. ( 3.11 )]. 

6*1° (t) and Ga-ir) are obtained directly from their expressions in terms of the fields f(j,e,p^ and d. As discussed 
in Sec. |l|, the chemical potential /i equals U/2 at half-filling and lies near the top of the LHB in the presence of a 
finite concentration of holes. The ground-state of the system in the atomic limit then has exactly one particle per site 
(even away from half-filling). We consider a single site and drop the site index. 

Let us first consider C'f^ (t) . It can be written as 

z(") 
1 



(^) ^ / '^^ I 2?[/,e,p,d]e*(r)pt(T)/T(r)/f(0)p|(0)e(0)exp{-5i^^} 

G'r(^) = -7M / '^^ 1 2?[/,e,P,d]d(T)p|(r)/^(r)/nO)pT(0)d*(0)exp{-5i"'^} 

dAzfgf(r), (CI) 



z(0) 
1 



where 



= j V[f,e,p,d\ exp{-5i°)}, 



9\°\r) ^ -— / P[/,e,p,d]e*(r)pTW/TW/f(0)pK0)e(0)exp{-5i?^}, 



Si M = / V[f,e,p,d]d{r)p^r)f,{T)fl{0)p^iO)d*{0)eM~Sl^'} (C2) 
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are the partition function and the Green's functions for a given configuration of the Lagrange multipliers, z^^^ = 
J dX z^"^ is the partition function for a single site. Since 5'^^^ is a Gaussian action when A is held fixed, we immediately 
deduce 

9f\r) = gf^iT)g,^i-T)g4T), (C3) 

where 

9fAr) = -e(^-^^^')-[^^(r)(l - F,) e{-r)F,], 

5e(r) = -e-^^*''- [0{t){1 + B) + 0(-r)i?] , 
g,,(r) = -e-(^^"'-<')-[^?(r)(l + i?T) + ^?(-r)i?T], 
9d{r) = -e-('^"'-^I^. ^' '+''^^ [0(r)(l + B') + e{-T)B'] , (C4) 
are the propagators of the fields fa, e, p] and d when the Lagrange multipliers are held fixed. We use the notation 

Fa = npi-H + iX'a^), 

B = nB{iX^^^), 

B' ^nsiiX^^^ -lY^Xl^^ +U), (C5) 



where np{x) — {e^^ + 1) ^ and nB{x) = {e^^ — 1) ^ are the Fermi and Bose factors, respectively. Using 



we obtain 

z(0) = / dX zf = 1 + e^^ + e/'C^M-c/) ^ (C7) 



in the limit T —t Q. In order to perform the integration over the Lagrange multipliers, we note that 

5^^g-»/3A(^)+./3A<^>^Q(g-2./3A<^))^ (C8) 

We deduce from Eqs. @, @ and (|C8|) 

where we have taken the limit T Q. The infinitesimal 77 (ry ^ 0+) is introduced to properly define the equal-time 
Green's functionsO 

Let us now consider the correction due to 

s.t-sl^^ ^ J dTA°j;p*f^p^. (CIO) 

Since Ga — z^^ J dXz\ga, the correction of order 0{A^) is given by 

Gi'' = + ^[ dXz^^^i^^ dXzfgi}\ (Cll) 



z^"' z^"' . z 
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where z^^^ and are the first-order corrections to zf^ and g'a^. The correction to the partition function is given by 



dX z 



(1) 



drAl 



From ( |CT^ ), we deduce 



The first-order corrections to gj'^'' and gj^"* are shown diagrammaticaUy in Figs. ^ and ||: 

dri A°(ti)5/^ (r - ri)g/^ (ri)gp^ (r - Ti)gp, (Ti)5e(-T) 
dri yl° (ti (r - n )gf^ (n )gp^ (0" )5p^ (T).ge (-r) 
+ y '^n^°(Ti)5/^(r)5p^(T-Ti)gp^(Ti).g/^(0")ge(-r), 
ff^W =y '^n^°(Ti)5/^(T)5p^(-T + Ti)g/^(0")5p^(-Ti).gd(T). 
Integrating over the Lagrange multipUers, we obtain 

-e^f^di-TMr - T,)e{~n) + e{-T + n)e{T,)]}, 

dAzf = e'^"""'" J dr, Al{T^){e^^^e{T)[e{T - T,)e{-T,) + e{-T + ti)0(ti)] 
From Eqs. (|cT^ ), (|Cll ) and ( |C16|) , we finaUy obtain 



dnGf\r-r,)A'^{n)G';'>in) 



(0), 



in the Umit T ^ 0. These results are easily extrapolated to higher orders as (in matrix form) 



(0) 



Gi = Gj 



GfAlGf 



Gfj^lGfAlGf 



or, equivalently, 



Gr^ = G 



(0)-l 

T 

(0)~1 



(C12) 
(C13) 

(C14) 



(C15) 



(C16) 



(C17) 



(C18) 



(C19) 



Gl^^Gf-^^Al. 

Note that these Green's functions, like the gauge field depend on the site which is considered. 
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FIG. 1. Diagrammatic representation of the contribution ^ 

of order t/U to the effective action S'lhb- (a) Quadratic 
part ^L^g. (b) Quartic part Slhb- When written in a 

time-ordered fashion, Sl^b generates a quadratic term which f* 



exactly cancels S'^^b (see text) 



(c) < ^ 



FIG. 2. Diagrammatic representation of the first-order 
{0{A°)) correction g\^^ to g\°> . The wavy line denotes the 
gauge field A°. 



fl 



> p^ 

> ^ 



d 



FIG. 3. Diagrammatic representation of the first-order 
(0(A°)) correction gf' to gf^. 
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